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Extensive efforts have been expended in devel¬ 
oping hybrid quantum systems to overcome the 
short coherence time of superconducting circuits 
by introducing the naturally long-lived spin de¬ 
gree of freedom. Among all the possible materi¬ 
als, single-crystal yttrium iron garnet has shown 
up very recently as a promising candidate for hy¬ 
brid systems, and various highly coherent inter¬ 
actions, including strong and even ultra-strong 
coupling, have been demonstrated. One dis¬ 
tinct advantage of these systems is that the spins 
are in the form of well-defined magnon modes, 
which allows flexible and precise tuning. Here 
we demonstrate that by dissipation engineering, a 
non-Markovian interaction dynamics between the 
magnon and the microwave cavity photon can be 
achieved. Such a process enables us to build a 
magnon gradient memory to store information in 
the magnon dark modes, which decouple from the 
microwave cavity and thus preserve a long life¬ 
time. Our findings provide a promising approach 
for developing long-lifetime, multimode quantum 
memories. 

Hybrid systems provide a promising solution for coher¬ 
ent information storage by combining the long coherence 
time of spin ensembles with the power of superconduct¬ 
ing circuits IDE]. Researches utilizing ensembles rang¬ 
ing from cold atomic gases [3] and magnetic molecules 
g to rare-earth-ion-doped crystals |SH7] and negatively- 
charged nitrogen vacancy (NV“) centres in diamond [8]- 
[12] have been reported recently. But not until very re¬ 
cently did people start to investigate the possibility of 
hybridizing yttrium iron garnet (YIG, Y 3 Fe 50 i 2 ), a ferri- 
magnetic insulator, with microwave cavities. In the YIG 
crystal, the coherent photon-spin ensemble interaction is 
greatly enhanced by the high spin density and can even 
approach the ultrastrong coupling regime p^HTS] . Gom- 
pared with previous dilute spin ensemble systems, direct 
spin-spin interaction leads to collective magnon excita¬ 
tions that have a range of distinct advantages, such as low 
damping rate, uniform distribution, rich nonlinear dy¬ 
namics, and well-defined mode profiles and wave-vectors. 
Particularly, each piece of YIG can be treated as a giant 
spin, with the flexibility to be individually manipulated 
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in experiments. 

In this work, we strongly couple multiple magnon 
modes with a microwave cavity resonance by placing mul¬ 
tiple YIG spheres into a three dimensional (3D) cav¬ 
ity. With these coherently coupled magnons, collective 
effects [16] such as a magnon dark mode (subradiant 
mode) that is decoupled from the environment, as well 
as enhanced interaction between a magnon bright mode 
(superradiant mode) and a microwave mode are demon¬ 
strated. The spectrum of the magnons can be further 
adjusted by controlling the local magnetic field of each 
YIG sphere, which allows us to engineer the dissipation 
and tailor the dynamics of the microwave photon at will. 
With this method, we applied a magnetic field gradient, 
which induced the periodic evolution of the magnons be¬ 
tween their temporal dark and bright modes, leading to a 
non-Markovian dynamics of the cavity energy that shows 
non-exponential decay and revival HD. By optimizing 
the coupling condition [18], good efficiency is obtained 
and our theory analysis indicates that a unity efficiency 
is achievable. Our experiment is performed at room 
temperature, demonstrating the coherent, long-lifetime, 
broadband and multimode gradient memory effect. Since 
the principle is entirely based on linear interactions, such 
a magnon memory can be readily scaled down to the 
quantum regime at millikelvin temperatures, providing 
a new approach besides the existing schemes for dilute 
spin ensembles [TM26) . This allows the realization of hy¬ 
brid quantum memories that do not suffer from common 
problems such as inhomogeneous broadening and so on. 

As a general situation, consider N identical YIG 
spheres loaded in a copper microwave cavity. The lin¬ 
ear coupling between the uniform magnon modes in the 
YIG spheres and the cavity TEno mode can be described 
by the Hamiltonian 

N 

PL/h = cjah'^'d -b + gjO^fhj + ^*dmj), (1) 

i=i 

where d, ifij are the bosonic operators associated with 
the photon and the uniform magnon modes, respectively, 
Qj is the coupling strength between the microwave mode 
and the magnon mode in the j-th YIG sphere with j 
ranging from 1 to A being the sphere index, uja is the 
resonance frequency of the TEno mode of the copper 
cavity with the YIG spheres loaded but at zero bias 
field, ujj = uJa P [j — {N P l)/2]Aco’^ are the evenly 
distributed magnon resonance frequencies, and is 
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their frequency interval. Such a frequency distribution 
is obtained by biasing each YIG sphere at a magnetic 
field Hj = i^o + [j -{N + l)/2]Ai7, where Hq = ujah 
is the common external magnetic field that brings a YIG 
sphere on resonance with the microwave cavity TEno 
mode, 7/27r = 2.8 MHz/Oe is the gyromagnetic ratio, 
and Ai7 = Acj^/y is the magnetic field difference be¬ 
tween neighboring YIG spheres that is provided by the 
fine tuning of the small coils underneath each individual 
YIG sphere (Fig.[^). 

We first study the coherence of a simple system with 
two YIG spheres (Y = 2). As a result of the mode hy¬ 
bridization for nonzero gi^2^ there exist three resonances 
in the cavity reflection spectrum when the two magnon 
modes are near resonance with the cavity TEno mode 
but not on resonance, which agrees well with our exper¬ 
imental observation (Eig.[^). By fixing AY and sweep¬ 
ing the external magnetic field around Yq, two avoided 
crossings are observed in the reflection spectra (Eig.[^), 
indicating the strong coupling between the cavity mode 
and the two magnon modes. When the frequencies of 
the two magnon modes (<^ 1 , 2 ) are brought closer by tun¬ 
ing AY, the absorption of the middle resonance becomes 
weaker, which eventually vanishes in the spectrum when 
the two magnon modes are brought simultaneously on 
resonance with the cavity mode {uji = UJ 2 = oja) by turn¬ 
ing off the gradient (AY = 0), and such a transition 
is evident in Eigs.[^-h. The remaining two resonances 
emerge from the hybridization of the bright (superradi¬ 
ant) magnon mode with the cavity photon mode, while 
the resonance that has disappeared from the spectrum 
is the dark (subradiant) magnon mode, as it decouples 
from the microwave cavity. 

The concepts of the bright and dark modes are il¬ 
lustrated schematically in Fig.[^. When the detuning 
field AY = 0, the two magnon modes are on resonance 
uji = UJ2- Since, in addition, we have gi = g2, the bright 
mode is the superposition of the two magnon modes that 
precess in phase B = ^ (mi + 7712 ), while the dark mode 
is the superposition of the two magnon modes that pre¬ 
cess out of phase D = ^ (mi —m 2 ). For the bright 
mode, the coherent interactions between the magnons 
and photons are collectively enhanced, leading to an en¬ 
hanced coupling strength g^ = V^gi, 2 ^ which is veri¬ 
fied in the avoided crossing spectrum where the splitting 
for AY = 0 {2gl27r = 18.88 MHz, Fig.[^) is times 
larger than that for AY = ±14 Oe (2^/27r = 13.42 MHz, 
Figs.[^&h). For the dark mode, the coupling of the two 
magnon modes with the cavity photons cancel each other, 
resulting in a vanishing coupling strength = 0- 

While the magnon bright mode can be used for in¬ 
formation conversion with microwave photons, the dark 
mode is an ideal candidate for information storage be¬ 
cause it decouples from the cavity and therefore has a 
very long lifetime m- A dark mode memory can be 
constructed if fast conversion between the bright and the 
dark mode is available, which in principle can be real¬ 
ized by rapidly tuning the magnetic bias field. However, 
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FIG. 1: Magnon dark mode of two YIG spheres, a, 

Device schematic including half of the copper cavity, two YIG 
spheres, and two small coils. Yq: external bias magnetic field; 
AY: magnetic field gradient generated by the small coils; h: 
magnetic field of the microwave resonance (TEno mode) in 
the copper cavity, b, Conceptual illustration of the magnon 
bright and dark modes: the bright mode is coupled to while 
the dark mode is isolated from the microwave magnetic field h. 
c, Reflection spectra as a function of AY under an external 
bias magnetic field Yo = 2858 Oe. d-h, Reflection spectra 
as a function of bias magnetic field Yo at a fixed magnetic 
field gradient AY = —14,-6,0,6,14 Oe, respectively. Each 
dashed line in c corresponds to the dashed line in each of the 
spectral maps in d-h. 


the slow response of the local inductive coils prohibits 
the experimental realization of such an approach. Al¬ 
ternatively, we consider here a magnon gradient memory 
(MGM) using magnon temporal dark modes which elim¬ 
inates the need of fast magnetic control. Such temporal 
dark modes are created by applying a magnetic field gra¬ 
dient (AY 7 ^ 0) to uniformly detune the magnon modes 
around the cavity resonance, so unlike the steady bright 
and dark modes discussed above, they are not the eigen- 
modes of the Hamiltonian. Instead, the non-Markovian 
dynamics of such a coupled system leads to a temporally 
evolving inter-conversion between the temporal bright 
and dark modes, which provides a good solution for fast 
memory operations without adding extra controls. 

For a device consisting of N YIG spheres, a frequency 
spectrum with Y ± 1 resonances spaced at an interval 
of Auj (~ Aujm for large Y) can be obtained as a result 
of the coupling between the cavity mode and the uni¬ 
formly distributed magnon modes. To achieve a memory 
with high efficiency, it is also important to ensure uniform 
coupling strengths for all the YIG spheres {gj = go, see 
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FIG. 2: The magnon gradient memory, a, Reflection spectra of the microwave cavity with eight YIG spheres at a bias 
held (Hq) of 2687 Oe (solid blue line) and 0 Oe (dashed black line), respectively. The frequency gradient of the magnon modes 
is tuned to be Au = 27r x 10 MHz. Inset: the rehection frequency spectra of the MGM at various bias magnetic he Ids. b, The 
dynamics of the MGM output for a 15 ns pulsed microwave excitation (gray) at the under-coupling (blue) and critical-coupling 
(red) conditions with a frequency gradient of Auj = 27r x 10 MHz. Dashed black curve shows the exponential decay when the 
magnons are largely detuned. Operation condition for the blue curve: u = 27rx 7.522 GHz, Hq = 2687 Oe; black: u = 27rx 7.522 
GHz, iTo = 0 Oe; red: cj = 27r x 7.120 GHz, Hq = 2544 Oe. The input frequency for each case is adjusted accordingly with 
the cavity resonance shift induced by the change of the external coupling. Inset: pulse retrieval at various frequencies for the 
under-coupling (top) and critical-coupling (bottom) condition, respectively. Zone I: rehection at the time of input due to the 
under-coupling condition of the external coupling; Zone H: hrst revival peaks; Zone HI: second revival peaks. 


Supplementary Information for details). The collective 
magnon modes from the N YIG spheres can be expressed 
as one temporal bright mode, expressed as 


7 = 1 


( 2 ) 


and N — 1 temporal dark modes, with the n-th dark mode 
expressed as 


1 ^ 

^ (3) 

V ^ j = l 

which do not couple to the cavity. Due to the 
uniformly-spaced detuning, the magnon modes described 
by Eqs. <§ and © will convert from one to another with 
a time interval T/N^ where T = 2tt jAuj is the evolu¬ 
tion period of the collective magnon modes. During the 
storing process, the cavity photons are hrst converted to 
the temporal bright magnon mode, which is then suc¬ 
cessively converted to the temporal dark magnon modes 
that prevent the radiation loss. After one period, the 
system evolves back to the bright mode, and the infor¬ 
mation retrieves from the magnons and converts back 
to microwave photons. Both storing and retrieving pro¬ 
cesses are accelerated due to the superradiance effect. In 
such a conhguration, no external control, such as refocus¬ 
ing pulse or gradient inversion, is required to retrieve the 


stored information. The performance of such an MGM 
depends on the number of YIG spheres. A larger N helps 
to improve the pulse re-construction and suppress the off- 
peak ripples (Supplementary Information). 

In our experiments, we demonstrate such an MGM us¬ 
ing eight YIG spheres {N = 8). The large tunability 
of the magnon allows us to engineer the cavity dissipa¬ 
tion with great flexibility and obtain a reflection spec¬ 
trum with a total of A" + 1 = 9 uniformly distributed hy¬ 
brid modes, as shown in Fig. Hf (solid blue line), where 
Aco’/27r = 10 MHz. Note that, in contrast to the sta¬ 
tionary dark mode, the temporal dark modes can be de¬ 
tected in the spectrum since they convert back to the 
bright mode periodically. In other words, each of the 
9 resonance lines in the steady state spectrum contains 
some component from the temporal bright mode which 
can be detected. When a 15-ns-long pulsed microwave 
signal at a frequency cj = = 7.52 GHz is injected into 

the cavity with an external bias magnetic field of 2687 
Oe, it couples to the magnon bright mode and then con¬ 
verts to the magnon dark modes. The retrieval of the 
stored pulse takes place after a pre-programmed time 
T = 21^1 Abj = 100 ns (blue peak in Fig.|^), in sharp 
contrast with the exponential decay when the magnons 
are strongly detuned by turning off the external bias 
magnetic field iTo (dashed black line in Fig.|^). While 
the MGM works without the requirement of any time- 
dependent control, extensions of the scheme would be 
possible when that control becomes available. For in- 
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stance, on-demand recall can be achieved via dynamic 
control by rapidly turning off and on the magnetic field 
gradient. Moreover, further dissipation engineering by 
controlling the field gradient would allow even more com¬ 
plicated manipulation of the output pulses, such as se¬ 
quence reversal, pulse splitting, and so on. However, a 
fast magnetic field tuning at this time scale with reason¬ 
ably large amplitude (around 50 Oe) is difficult to achieve 
and so it falls beyond the scope of this work and is left 
for future study. 

The performance of the MGM can be significantly im¬ 
proved by optimizing the external coupling condition. 
Blue curves in Figs. & b correspond to the condition of 
under-coupling and therefore the retrieved pulse is very 
weak. A drastic boost of the revival pulse is obtained 
(red peak in Fig. It) by adjusting the coaxial probe to 
meet the critical coupling (impedance matching) condi¬ 
tion: Ka,i = /^a,o+'7r|^p/Aco’, where Ka,o is the amplitude 
damping rate of the cavity resonance and Ka,i is the cou¬ 
pling rate from the coaxial probe. The critical coupling 
condition allows a complete energy conversion into the 
magnon bright mode (see Supplementary Information), 
and it also results in a maximized signal retrieval back 
to the microwave photon. Therefore, the first revival 
peak will contain all the stored information and there will 
be no successive revival peaks as observed in the under- 
coupling situation. Such analysis matches our measure¬ 
ment results shown in the insets of Fig. It: the second re¬ 
vival peaks (zone III) are relatively strong compared with 
the first revival peaks (zone II) for the under-coupling 
situation (top inset), while for the critical-coupling situ¬ 
ation (bottom inset) the second revival peaks disappear. 

One advantage of the MGM is that its operation is not 
restricted to a narrow frequency range. A uniform fre¬ 
quency distribution exists in a wide frequency range as we 
vary the bias magnetic field (Fig.[^, inset). Accordingly, 
broad-band pulse retrieval process can be obtained in 
the time-domain measurement. Generally, the operation 
bandwidth can be determined by NAuj, and therefore, 
the delay-bandwidth product is set by N. In the exper¬ 
iment, we observed retrieval peaks in a broad frequency 
band of 80 MHz at a bias magnetic field of 2687 Oe 
(Fig-IJj, insets), which is 1-2 orders of magnitude larger 
than the linewidth of the photon or magnon resonances. 
Note that when sweeping the input frequency, the revival 
peak intensity periodically varies for the under-coupling 
situation (Fig. It, top inset, zone H), which can be at¬ 
tributed to the periodic variation in the coupling condi¬ 
tion, while such a phenomenon disappears in the criti¬ 
cally coupled situation thanks to the improved coupling 
condition. 

We quantitatively characterized the MGM by measur¬ 
ing the delay time and the retrieval efficiency as a func¬ 
tion of Acj, and the result is plotted in Fig.[^. As Auj re¬ 
duces, the revival time T increases, and the measurement 
results perfectly follow the relation T = 27r/Aco’ (Fig.[^). 
The MGM retrieval efficiency is defined as the ratio of 
the output pulse energy to the input pulse energy. It has 
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FIG. 3: Characterization of magnon gradient memory. 

a, Measured retrieval pulses for different frequency gradients 
Aa;/27r. The input pulse is only shown for Act;/27r = 10 MHz 
and remains the same for the other cases, b, Extracted stor¬ 
age time as a function of frequency gradient. The solid line 
is calculated by T = 27 r/Act;. c, Retrieval efficiency as a 
function of the storage time obtained from the measurement 
(circles), the numerical fitting (solid blue line), and the calcu¬ 
lation using Eq. 0 (dashed black line), respectively, d, Mea¬ 
sured interference between the retrieved and reference signals 
as a function of their phase difference. The solid line is a 
sinusoidal fitting. 


been extracted as a function of delay time and plotted in 
Fig.§ (blue circles), showing an efficiency of about 30% 
under the critical coupling condition. The extracted effi¬ 
ciencies match the numerical simulation results obtained 
using the same parameters as in the measurement (solid 
blue line). On the other hand, the efficiency can be ap¬ 
proximated using the asymptotic expression (see Supple¬ 
mentary Information) 

G 

TTg 

which is plotted in Fig. 0 (dashed black) for compari¬ 
son. Equation @ shows that the MGM efficiency de¬ 
pends on the figure of merit G = | ^ for an input pulse 
with a duration of where F = is the finesse of 

the magnon gradient and G = ^ ^ is the magnon- 
photon interaction cooperativity. Here we assume all the 
magnon modes in each YIG sphere have the same am¬ 
plitude damping rate (/^m)- Therefore by optimizing G 
and F we can further improve the MGM performance. 
In practice, F can be well controlled to be in the range 
of 10^ — 10^, and G in the order of 10^ — 10^, and thus 
C ~ 1 is achievable for real applications. 

In addition to the high achievable efficiency, such an 
MGM preserves the signal coherence during the storage 
and retrieval process. Figurej^ plots the measured inter- 
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FIG. 4: Multi-pulse storage in the MGM. Retrieved 
pulses for a double-pulse excitation. The two pulses are sepa¬ 
rated by 40 ns with a 15 ns duration each. Inset: double-pulse 
retrieval for various input frequencies at a bias held of 2687 
Oe. 

ference signal between the retrieved pulse (after a 100-ns 
storage time) and the input microwave signal. The peri¬ 
odic dependence of the interferometer output on the rel¬ 
ative phase between these two signals indicates they are 
still phase coherent. A visibility as high as 96.9 ± 1.5% 
can be extracted, proving the high coherence preserva¬ 
tion feature of our MGM. 

Another advantage of the MGM is its multimode op¬ 
eration capability, which is crucial for speeding up quan¬ 
tum computation protocols. Figure]^ plots the output 
of the MGM when two identical pulses are sent into the 
memory. The input pulses are 15 ns in duration and sep¬ 
arated by 40 ns. Two pulses are retrieved from the mem¬ 
ory after the pre-programmed storage time T = 100 ns. 
As there is no echo process involved during the memory 
process, the MGM operates in a first-in-first-out (FIFO) 
mode. The multimode operation can be also achieved 
in a wide frequency range (Fig.[^ inset). The occasional 
pulse distortion at certain frequencies is a result of the 
non-ideality of the system such as the non-identical cou¬ 
pling strengths and decay rates of different YIG spheres, 
or the impact of the high order magnon modes. 

In summary, we demonstrated coherent coupling of 
multimode magnon resonances to a single microwave cav¬ 
ity. The dark magnon mode of two YIG spheres is ob¬ 
served. We used the temporal magnon dark mode to es¬ 
tablish non-Markovian dynamics in a system consisting 
of eight YIG spheres, and showed its potential in devel¬ 
oping quantum memories, which are broadband, multi- 
mode, and do not require fast field switching. Our inves¬ 
tigation on a novel mechanism of manipulating magnons 


and microwave photons paves the route towards magnon- 
based hybrid quantum memories, quantum repeaters as 
well as quantum networks [28] . 

Methods 

Sample preparation. The microwave cavities are ma¬ 
chined from high-conductivity copper. The cavity hous¬ 
ing two YIG spheres has a dimension of 43 x 21 x 9.6 mm^, 
with its TEiio resonating at 7.87 GHz, while the one 
housing eight YIG spheres has a dimension of 50 x 21 x 4.8 
mm^, with its TEno mode resonating at 7.52 GHz. All 
the YIG spheres are highly polished and have a diameter 
of 0.25 mm, with their (110) direction aligned parallel to 
the external bias magnetic field. The mounting ceramic 
rods that hold the YIG spheres are glued on the cavity 
wall, close to the the maximum magnetic field of the cav¬ 
ity resonance to achieve strong magnon-photon coupling. 
Eight holes, 6.5 mm apart from each other, are drilled 
deep into (but not through) the cavity wall underneath 
the YIG spheres, where eight home-made small coils (200 
turns, 6 mm in diameter, made of 32-gauge copper wires) 
are inserted with a 1 mm distance from the YIG spheres 
for fine tuning of the bias magnetic fields within a range 
of ±30 Oe (Eigs. la & S5a). The whole device assembly is 
inserted between the two poles of a water-cooled electro¬ 
magnet which provides the strong external bias magnetic 
field. 

Microwave characterization. The frequency spectra 
are taken using a vector network analyzer by measuring 
the reflection signal from the coaxial probe which ac¬ 
cesses the cavity mode through a small hole in the cav¬ 
ity wall. The external coupling can be tuned by adjust¬ 
ing the position of the coaxial probe. The time traces 
are measured using a high-speed oscilloscope when the 
continuous-wave (CW) input microwave signal is modu¬ 
lated by a pulse generator through a transistor-transistor 
logic (TTL) switch which gives a rectangular pulse shape. 
To measure the coherence between the retrieved pulse 
and the input signal, the CW input signal is split into 
two branches before the pulse modulation, one of which, 
serving as the reference, is phase adjusted and combined 
with the retrieved signal before sending into an envelope 
detector. Then, the visibility of interference between the 
retrieved signal and the reference is measured by varying 
the reference phase. More details about the microwave 
characterization can be found in the Supplementary In¬ 
formation. 

Theoretical derivation. In the Supplementary Infor¬ 
mation, details about the theory derivation and analysis 
of the dark mode and the gradient memory are provided. 


neering. Phys, Scr. T137, 014001 (2009). 

[2] Xiang, Z.-L., Ashhab, S., You, J. & Nori, F. Hybrid 


[1] Wallquist, M., Hammerer, K., Rabl, P., Lukin, M. & 
Zoller, P. Hybrid quantum devices and quantum engi- 





















6 


quantum circuits: superconducting circuits interacting 
with other quantum systems. Rev. Mod. Phys. 85, 623- 
653 (2013). 

[3] Verdii, J. et al. Strong magnetic coupling of an ultracold 
gas to a superconducting waveguide cavity. Phys. Rev. 
Lett. 103, 043603 (2009). 

[4] Eddins, A. W., Beedle, C. C., Hendrickson, D. N. & 
Friedman, J. R. Collective coupling of a macroscopic 
number of single-molecule magnets with a microwave cav¬ 
ity mode. Phys. Rev. Lett. 112 , 120501 (2014). 

[5] Schuster, D. I. et al. High-cooperativity coupling 
of electron-spin ensembles to superconducting cavities. 
Phys. Rev. Lett. 105, 140501 (2010). 

[6] Probst, S. et al. Anisotropic rare-earth spin ensemble 
strongly coupled to a superconducting resonator. Phys. 
Rev. Lett. 110, 157001 (2013). 

[7] Tkalcec, A. et al. Strong coupling of an Er^^-doped 
YAIO3 crystal to a superconducting resonator. Phys. 
Rev. B 90, 075112 (2014). 

[8] Zhu, X. et al. Coherent coupling of a superconducting 
flux qubit to an electron spin ensemble in diamond. Na¬ 
ture 478, 221-224 (2011). 

[9] Kubo, Y. et al. Hybrid quantum circuit with a super¬ 
conducting qubit coupled to a spin ensemble. Phys. Rev. 
Lett. 107, 220501 (2011). 

[10] Amsiiss, R. et al. Cavity QED with magnetically cou¬ 
pled collective spin states. Phys. Rev. Lett. 107, 060502 
( 2011 ). 

[11] Marcos, D. et al. Coupling nitrogen-vacancy centers in 
diamond to superconducting flux qubits. Phys. Rev. Lett. 
105, 210501 (2010). 

[12] Ranjan, V. et al. Probing dynamics of an electron-spin 
ensemble via a superconducting resonator. Phys. Rev. 
Lett. 110, 067004 (2013). 

[13] Huebl, H. et al. High cooperativity in coupled microwave 
resonator ferrimagnetic insulator hybrids. Phys. Rev. 
Lett. Ill, 127003 (2013). 

[14] Zhang, X., Zou, C.-L., Jiang, L. & Tang, H. X. Strongly 
coupled magnons and cavity microwave photons. Phys. 
Rev. Lett. 113, 156401 (2014). 

[15] Tabuchi, Y. et al. Hybridizing ferromagnetic magnons 
and microwave photons in the quantum limit. Phys. Rev. 
Lett. 113, 083603 (2014). 

[16] Dicke, R. Coherence in spontaneous radiation processes. 
Phys. Rev. 93, 99-110 (1954). 

[17] Zou, C.-L. et al. Photonic simulation of system- 
environment interaction: non-Markovian processes and 
dynamical decoupling. Phys. Rev. A 88, 063806 (2013). 

[18] Afzelius, M. & Simon, C. Impedance-matched cavity 
quantum memory. Phys. Rev. A 82, 022310 (2010). 

[19] Hau, L. V., Harris, S. E., Dutton, Z. & Behroozi, C. H. 
Light speed reduction to 17 metres per second in an ul¬ 
tracold atomic gas. Nature 397, 594-598 (1999). 

[20] Anderson, A. et al. Spin echo serial storage memory. J. 
Appl. Phys. 26, 1324-1338 (1955). 

[21] Alexander, A. L., Longdell, J. J., Sellars, M. J. & Man- 
son, N. B. Photon echoes produced by switching electric 
fields. Phys. Rev. Lett. 96, 043602 (2006). 

[22] Hosseini, M., Sparkes, B. M., Campbell, G., Lam, P. K. 
& Buchler, B. C. High efficiency coherent optical mem¬ 
ory with warm rubidium vapour. Nat. Commun. 2 , 174 
( 2011 ). 

[23] Hedges, M. R, Longdell, J. J., Li, Y. & Sellars, M. J. 
Efficient quantum memory for light. Nature 465, 1052- 


1056 (2010). 

[24] Kraus, B. et al. Quantum memory for nonstationary 
light fields based on controlled reversible inhomogeneous 
broadening. Phys. Rev. A 73, 020302 (2006). 

[25] Arnold, S., Whitten, W. B., Ramsey, J. M. & Liu, C. T. 
Room-temperature microparticle-based persistent spec¬ 
tral hole burning memory. Opt. Lett. 16, 420-422 (1991). 

[26] de Riedmatten, H., Afzelius, M., Staudt, M. U., Simon, 
C. & Gisin, N. A solid-state light-matter interface at the 
single-photon level. Nature 456, 773-777 (2008). 

[27] Dong, C., Fiore, V., Kuzyk, M. C. & Wang, H. Optome¬ 
chanical dark mode. Seienee 338, 1609-1613 (2012). 

[28] Kimble, H. J. The quantum internet. Nature 453, 1023- 
1030 (2008). 

Acknowledgments 

This work is supported by DARPA/MTO MESO pro¬ 
gram. H.X.T. acknowledges support from a Packard 
Fellowship in Science and Engineering. L.J. acknowl¬ 
edges support from the Alfred P. Sloan Foundation, 
the Packard Foundation, the Multidisciplinary Univer¬ 
sity Research Initiative (MURI), and the DARPA Qui- 
ness program. C.L.Z. acknowledges National Basic Re¬ 
search Program of China (Grant Nos. 2011CB921200 
and 2011CBA00200). F.M. acknowledges ERC OP- 
TOMECH and DARPA ORCHID. The authors thank 
Dr. Michel H. Devoret and Dr. Michael Hatridge for 
providing a prototype 3D microwave cavity. 

Author contributions 

X.Z. prepared the samples, performed the measurements 
and data analysis. C.L.Z. and L.J. conceived the idea. 
C.L.Z., L. J. and F.M. provided theoretical analysis. N.Z. 
assisted with the measurements. L.J. and H.X.T. super¬ 
vised the project. X.Z., C.L.Z., F.M., L.J., and H.X.T. 
wrote the manuscript. 

Additional information 

Supplementary Information is available in the 
online version of the paper. Reprints and 

permissions information is available online at 
www.nature.com/reprints.Correspondence and re¬ 
quests for materials should be addressed to H.X.T. 
(hong.tang(Ryale.edu). 

Competing financial interests 

The authors declare no competing financial interests. 



so 

Supplementary Information 

Contents 

References [5] 

I. Microwave photon-magnon strong coupling [ST] 

II. Reflection spectrum [ST] 

III. Dark magnon modes [S2] 

A. Two Magnon Modes [S2] 

B. Generalization to Multiple Magnon Modes [S2] 

C. Temporal Dark Magnon Mode [S3] 

IV. Magnon gradient memory [S3] 

A. Intuitive Explanation [S3] 

B. Dynamics [S4] 

C. Asymptotic Solution [S4] 

D. Scattering Picture [S5] 

E. Critical coupling condition [S6] 

E. Photon Retrieval [S7| 

G. Measurement Scheme [S8] 

H. Experimental Imperfections [S8] 


SI 


I. MICROWAVE PHOTON-MAGNON STRONG COUPLING 


The interaction between microwave photon and magnon (single mode) can be described by the Hamiltonian with 
the rotating-wave approximation 

+ ^(a^m + am^), (SI) 

where d) (a) is the creation (annihilation) operator for the microwave photon at frequency cJa.For the uniform magnon 
mode, the collective spin excitations are approximately represented by the Boson operator m^{m) with Holstein- 
Primakoff approximation [SI]and ujrn is the magnon frequency. The coupling strength g between the two systems is 
[S2]: 




V2Ns, 


(S2) 


where uj is the resonance frequency, 7 is the gyromagnetic ratio, Va is the modal volume of the microwave cavity 
resonance, go is the vacuum permeability, N is the total number of spins, and 5 = | is the spin number of the ground 
state Fe^+ ion in YIG. The coefficient 77 < 1 describes the spatial overlap and polarization matching conditions 
between the microwave field and the magnon mode, which can be explicitly written as: 

(t(r)-7)" + (t(r).7)2 

77 — - ^ -, (^007 

max{| h (r)p} 


where max{| h (r)p} is the maximum magnetic field intensity of the cavity mode, and h (r) is the magnetic field 
amplitude at the location (r) of the YIG sphere, ^ with j = x^y^z are unit vectors and ^ is along the bias field 
direction. 


II. REFLECTION SPECTRUM 


In our experiment, we measure the reflection spectrum under continues driving and the emission after pulsed 
excitations of the microwave cavity, which correspond to the eigenmode spectrum and the transient cavity ouput, 
respectively. In general, the equation of motion for the cavity photon is 

= -i[a, H] - (S4) 

where l-L is the system Hamiltonian, tZag is the cavity external coupling rate to the coaxial probe, Ein{t) and uji are 
the amplitude and frequency of the input microwave. 

In the rotating frame of cj/, the microwave electric field collected by the detector reads 

Kut(i) = -Ein{t) + i^2Ka,ia{t) (S5) 

For a continuous stable input Ein{t) = Ein, we measure the stationary spectra with the normalized spectrum 
expressed as 

r{uji) = -1 + (S6) 

L/in 

For the memory operation, the retrieval of the microwave pulse is measured after the input pulse is stored, and the 
detected intensity is 

I{t)=2Ka,l\m)\^- (S7) 


The efficiency of the MGM, defined as the energy ratio of the retrieval pulse to the input pulse, can be expressed as 


^ f I{t)dt 

J Iin{t)dt' 


(S8) 


[51] Holstein, T. and Primakoff, H. Phys. Rev. 58, 1098 (1940). 

[52] Zhang, X., Zou, C.-L., Jiang, L., and Tang, H. X. Phys. Rev. Lett. 113, 156401 (2014). 



S2 


III. DARK MAGNON MODES 

In general, the Hamiltonian of multiple linearly coupled magnon modes and the cavity mode reads: 

N 

'H/h = co’ad'^'d + + gjd)mj + (S9) 

i=i 

When deriving the equations of motion, both the magnon dissipation and the cavity dissipation (including the external 
coupling) are considered. In principle, the input-output theory shows that this gives a damping term as well as a 
fluctuation (including vacuum and thermal) term. To keep the notation simple, we will, from now on, use symbols 
like a and rrij to refer to the expeetation values of the underlying operators. Therefore, the fluctuation terms will be 
dropped. This is sufficient to obtain the linear scattering behavior of the system and to discuss its eigenmodes. With 
this proviso, the equations of motion for the cavity mode and the magnon modes are: 

fa = {-icoa-Ka)a-i^gjmj 

3 

frrij = {-iwj - Kj)mj - ig*a. 

Here, Ka and kj are the total decay rates of the cavity photon and the magnon modes, respectively, with Ka including 
both the internal dissipation (/^a,o) as well as the coupling to the coaxial probe (/^a,i)- 


(510) 

(511) 


A. Two Magnon Modes 

As described in the main text, we first study the case with two magnon modes. For simplicity, let the two magnon 
modes be on-resonance with the cavity mode (cji = 002 = ^m)- Then the dynamics of the cavity field is 

fa = {-ioJa - Ka)a - i{gimi + g2m2). (S12) 

The cavity mode couples to the collective mode B = {givai + g2m2)/\J\gi |V| 52 |'as 

fs = {-iuJm - Km)B - iflgil'^ + |52pa. (S13) 

And the orthogonal mode D = {g^mi — gim2)/f \gif‘ + \g2\t‘ satisfies 

fo = i-iuJm - Km)D - i{glg2 - 929*i)a/\j\gi\^ + \g2\^ (S14) 

— ( iojjY^ Kjjyi)D. (S15) 


Therefore, D is completely isolated from the cavity mode, representing the dark mode that cannot be detected from 
the reflection spectrum of the microwave cavity. 


B. Generalization to Multiple Magnon Modes 


For multiple magnon modes N >2 and = cji for j = 2,..., A, the normalized bright mode is 

N 

B 

i=i 

Then we have 




(SI 6 ) 


(S17) 




If the magnons are identical, the bright mode is a collective mode of all the magnon modes with the coupling strength 
being enhanced by a factor of 


/ = = x/iv. 

There are N — 1 other modes, orthogonal to the bright mode and decoupled from the cavity. 


(S18) 


C. Temporal Dark Magnon Mode 

In the above analysis, all the magnon modes are on resonance with the cavity mode, and therefore the bright and 
dark modes are all the eigenmodes of the Hamiltonian. When the magnon modes are tuned off but close to resonance, 
situation would be different. As an example, we consider the case of two magnon modes. If the two magnon modes 
are detuned from the cavity mode as 


uJi = uJa^ Aa;/2, 

UJ2 = uJa - Aa;/2, 

then we can define the temporal bright mode at time t as (assuming gi = g2 = Qm) 


B«) = 


‘'at 




■ 


(519) 

(520) 


(S21) 


At t = 0, this mode is B = + m 2 ), which is the same as static bright mode discussed above. But at t = 

the mode evolves to 




(S22) 


which is the dark mode that decouples from the cavity. Therefore, in this scenario, the temporal bright and dark 
magnon modes are not eigenmodes of the system and they are inter-convertible through time evolution. 


IV. MAGNON GRADIENT MEMORY 


A. Intuitive Explanation 


To realize a magnon memory device using the temporal magnon dark mode, we want: (1) the input photon can 
be converted to the collective magnon bright mode as quickly as possible, and therefore the photon energy will not 
dissipate too much due to the cavity intrinsic loss; (2) the bright mode can convert to the dark mode once the input 
photon is converted to the bright mode; (3) the retrieval of the photon is predictable and pre-programmable. 

To fulfill these requirements, the system should consist of multiple magnon modes with equal detuning. Assuming 
the lifetime of the magnon and the cavity photon are = 1/2/^^ and = 1/2/^^, respectively, the retrieval period 
T = tt/Acj, the conversion time from the bright mode to the dark modes = T/N with N being the total number 

of the magnon modes, and the conversion time from the cavity photon to the magnon bright mode ta^B = 
then we will have the following constraints 


Acj 

9 


T < Tm, 

(S23) 


(S24) 

ta^B ^ "Ta? 

(S25) 

TT 

> - = 27r X 

(S26) 



> ^/N/S.UJ > 2Tl\fN X Km, 

(S27) 


27r 


which can be rewritten as 
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A highly efficient magnon gradient memory requires the cooperativity C = —-— ^ 47r^. In our experiment, the 

hybrid YIG sphere-3D cavity structure can have a cooperativity C as large as 10^, and therefore it is very promising 
for memory applications. 


B. Dynamics 


In an ideal magnon gradient memory, we have 


N -1 


A(j + jAuj 


for the uniform magnon mode in the j—th. YIG sphere. Since 


^rrij = {-iwj - Kj)mj - ig*a, 


we have the formal solution for the magnon time-evolution (where mj{0) = 0): 


mi{t) = -ig* [ 

Jo 


Then, the cavity dynamics is 


d 

— a 
dt 


= {-icoa-Ka)a-y^\g/ 

)a- [ a(T)y'|gjf 

Jo 


= (-■- Ka)a - 


For N identical YIG spheres with gj = g and K,j = Km, we have 

3 3 


N 


— |^|2 ^{-iuJa-i^^^^Auj-Km){t-r) ^-ijAuj{t-T) 




Finally, the cavity dynamics is determined by the following equation 


sm 


j=i 

N Auj{t—r) 1 


sm 


Au;(t—r)-| 


d 

■^a{t) = {-ioJa - K.a)a{t) - \gf' J a(T)( 




sm 


NAuj{t—r) -I 


sm 


Auj{t—r) - 


-dr. 


(529) 

(530) 

(531) 

(532) 

(533) 


(S34) 


(S35) 


C. Asymptotic Solution 


As shown by Fig.|Sl[ we can take the approximation 


sm 


NAuj{t—r) -I 


3^ ss 27ry](-l)'(^ ^^6[Aoj{t-T -l^)] 


sm 


Auj{t — r) 


Auj ^ 


for Y ^ 1, where 5{x) is the Dirac delta function. 

Therefore, we can describe the cavity dynamics for Y ::$> 1 as 


d 

dt 


(S36) 


pt ^ o 

a{t) = [-iuja-Ka)a{t)-2'K\gf‘ j - t - l — )]dT 

^ l= — oo 


•i^r 




= ( — iuJa — Ka)a(t) -, , 

V a a! K I Aw ^ ' Aw ' Aw' 

Ir^O 


(S37) 
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FIG. SI: The plot of the function f{x) = 

Sin [ ^ Xj 


For t < the cavity photon experiences single exponential decay 


d , ^ , X / X tt IqI 

= {-luja - na)ci{t) - 


with solution 


a{t) = a(0)e 




We can see that the cavity photon decay is enhanced by 


For i < 2|^, we have 


d ,, 27r|(7p , 27r , 

-ait) = + a(i-—)e(- - - 


= (-tWa - Ka 


•i-ir 


)a{t) + 


271 i^r 


a(0)e^ 


( - iWa - Km ) g-(iWa+fCa + 7=^ ) (t-S5r) 


Aw Aw 

Therefore, the cavity photons show a revival due to the energy of magnons coupling back to the cavity 


ait) = - 




(S38) 


(S39) 


(540) 

(541) 

(542) 


D. Scattering Picture 


Above, we have adopted a direct solution of the temporal equations of motion. Alternatively, one can also d iscus s 
the MGM as a scattering problem in frequency space. We employ the equations of motion, Eqs. (SIO) and (Sll), 
supplement them with the input field (as in Eq. ([S4|)), and obtain, in frequency space: 


(543) 

(544) 


-iuja = {-iuja - i^a)ci - i ^ gjrrij - 

j 

—iujrrij = {—iuJj — Km)Tnj — ig^a . 

Here the functions a, m^, and E[^ are understood to be functions of the frequency uj. 

These equations can be solved by eliminating mj in favor of a and then applying the input-output relation Eq. (S5). 
This leads to the following result for the frequency-dependent reflection amplitude that gives £’out(<^) = 
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r(w) = 


[Wa-U + S(a;)] + i{Ka,i - Kgfi) 

— [Wq — U! + S(w)] + i{Ka,l + Ka,o) 


(S45) 


At this point we have introduced the “MGM self-energy” that describes the collective effects of all the magnon modes 
acting on the cavity mode: 


= E 

3 


UJ — UJj 


(S46) 


The real part of S describes the effective frequency shift, while the imaginary part describes additional damping 
(induced by the magnon mode damping). The reflection amplitude can be decomposed into magnitude and phase 
shift: 


rH = |r(w)|e*'’(‘"). (S47) 

In the ideal case without intrinsic losses (/^a,o = = 0), we have \r\ = 1. The time-delay of a scattered wave-packet 

is determined by the derivative of the phase shift with respect to frequency: 


= (S48) 

The most important necessary condition for a useful memory is that this time-delay be constant over the bandwidth 
interval NAuj, i.e. the slope of 0 should be constant. In the limit of small magnon-cavity coupling ^ ^ 0 (and no 
intrinsic losses), each magnon resonance leads to a step of 27r in the phase shift 0. A finite g rounds off these steps. 
In the vicinity of each resonance, we can then approximate T, ^ g‘^/{uj — ujj)^ which leads to a phase shift rising like 
0{uj) — 0{ujj) = 2i<ia^i{o3 — ujj)/g‘^. In order to have a constant overall slope of 6>, we have to match this to the slope 
T = 27r/Aco’ that is dictated by the spacing of resonances and which corresponds to the ideal storage time. That leads 
to the critical coupling condition (without intrinsic losses) 


Auj 


= . 


(S49) 


We note that this condition remains true (in the present form) if intrinsic cavity losses are also incorporated, i.e. 
when ^ 0. The contribution from /^a,o cancels when deriving the condition. We note, however, that for finite 
/^a,o this condition slightly differs from the critical coupling condition derived below from demanding zero reflection 
of the input pulse. This is because, in general, for finite /^a ,0 5 the form of the critical coupling condition depends on 
the precise physical condition that is imposed. Still, one needs ideally /^a,i to suppress unwanted losses and 

to avoid a resonance structure showing up in the magnitude |r| of the reflection. In addition, the cavity mode acts 
like a filter, which should be broad enough to cover the whole magnon spectrum, i.e. 


^ NAuj . (S50) 

Taking these two conditions together also implies g /^a,i/V^, as stated already in Eq. ( |S28|). In addition, to keep 
the magnon losses small (ensuring |r| close to 1), one needs Hm <C Auj j21^^ as stated in Eq. (|S26|. 


E. Critical coupling condition 


Erom Eq. (S39), the microwave pulse input into the cavity couples to the external coupling channel (rate Ka,i) 
and the temporal bright mode ( ^^^^ ), and also dissipates (/^a,o) due to intrinsic radiation and absorption losses. 


Therefore, the full internal loss of the cavity is actually K,a o ~ ^a,o 


Au; 


To suppress the reflection of the input 


pulse, the critical coupling (impedance matching) condition = ^a,o ~ desired. The storage and 

retrieval efficiency of the MGM can be greatly affected by the critical coupling condition of the MGM with the coaxial 









S7 





FIG. S2: Input reflection at various coupling conditions, a, Measured input reflections for microwave pulses to the microwave 
cavity with under-coupled (ACa,i < ^a,o? blue curve) and critically coupled (/€a,i ~ ^a ,0 5 I'Gd curve) probes, respectively, b, The 
calculated cavity reflection for a rectangular pulse input (shadow region) at various coupling conditions, c, The calculated 
cavity reflection for a Gaussian input pulse using the same parameters as in b. 


probe. If the MGM is under-coupled, a large portion of the energy will be reflected as the signal gets in and out of 
the cavity, and as a result the MGM efficiency becomes low. 

In our measurements, we did observe a strong reflection peak at the time of input for the undercoupling situation, 
and such reflection peak can be eliminated after adjusting the external coupling rate to meet the critical coupling 
condition, as plotted in Fig. S2r. Such an observation agrees well with the theoretical predication (Fig.[S^). Note 
that there is still some residual reflection for the critical coupling condition, which is attributed to the step response 
of the MGM to the rectangular input pulse (both the rising and falling edges). This imperfection can be suppressed 
by choosing Gaussian pulse inputs (Fig.IS^). 


F. Photon Retrieval 


From Eq. (S42), after one period T = ^ the stored energy will couple back to the cavity, which leads to the retrieval 
of the photons. We can solve for the detected microwave amplitude by the input-output formula Eq. (S5), with the 
input microwave also taken into account. Substituting the experiment parameters into the equation of the cavity 


field, with the input rectangular pulse duration tp = 20 ns, we obtained the dynamics in Eig. S3r. The analytical and 
numerical results agree well with each other. With the given parameters, the numerical solution shows a memory 
efficiency of C = 0.33, which is the integrated total output energy in the first retrieval peak as compared with that of 
the input pulse. 

Erom the asymptotic solution, we can derive the efficiency of the magnon gradient memory as 


— ‘2k 

g Auj 


1 - 


tp ) \^a,0 + 


zeM! 

Au! 


7 ^\ 9 \ 

Au! 


(S51) 


for a rectangular input pulse under critical coupling condition (which demands zero reflection). Note that the pulse 

duration is limited by the external coupling rate: tp > = l/(/^a,o + 

In Eig.[S3|3, we plot the efficiency as a function of the coupling strength go. A comparison of the numerical results 
with the asymptotic solutions shows that they agree well with each other for 5f/27r > 7 MHz. The efficiency saturates 


to C = e 


0.403 for 


^\9oV 

Au) 


^ which is limited by the intrinsic loss of the magnon. Eigure S3" plots the 


behavior of the memory for various coupling strengths, which shows the saturated retrieval peak at 100 ns. If we can 
reduce the magnon dissipation rate, as shown in Eig.[S3|i, the retrieval efficiency approaches unity. If operating at 
low temperatures, the magnon linewidth can be reduced to 0.042 MHz [S3], which will lead to a saturated efficiency 
of 0.95 for a storage time ^ = 100 ns. 


[S3] Spencer, E. G., LeCraw, R. C., and Linares, R. C., Jr. Phys. Rev. 123, 1937 (1961) 
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By numerical calculation, we further studied the memory for different input pulse detuning and frequency gradient 
of the magnons. From our analytical solutions, we can expect the effective bandwidth of the MGM to be about 
^ao — ^a,o + ~ for ^ oo. Howcvcr, in experiments we have a finite and therefore the bandwidth is 

also limited by the bandwidth of the gradient magnon spectrum, which is about NAuj. Therefore, the highly efficient 
memory works in a bandwidth of min{ ^|^^^ , NAuj}. From Fig.|S4^, the bandwidth is about SAcj, which is consistent 
with our theory. Figure[S4)3 shows that the storage time is inversely proportional to Acj, which also agrees with our 
expectation. 


G. Measurement Scheme 


FigurejS^ illustrates the device schematic of the MGM with eight YIG spheres. Note that the radius (3 mm) of the 
small coils is much larger than the YIG sphere radius to ensure the magnetic fields generated by the coils are uniform 
at the position of the YIG sphere (1 mm above the coil). As a result, the coils have to be placed 6.5 mm apart from 
one another, and so are the YIG spheres. Since the magnetic field of the cavity TEno mode has a cosine distribution 
along the direction inside the cavity (Figs. S5 3 & c), different YIG spheres experience different magnetic field 
strengths. Therefore the coupling strength of the magnon modes with the cavity mode are not identical, and this 
contributes to the non-ideality of the MGM. Such non-ideality is unavoidable, but in our experiment efforts such as 
placing the YIG spheres as close to the center as possible have been taken to reduce the non-ideality in the device. 

There are three types of measurement involved in our experiments. 

(1) First is the cavity reflection spectrum measurement. The input signal is provided and the reflected signal is 
detected by a vector network analyzer, as indicated in Fig.[S6^. A circulator is used to separate the input and reflected 
signals to avoid undesired interference. 

(2) The second type of measurement is the time trace measurement^ which is carried out using a high-speed oscillo¬ 
scope (Fig.|S6|3). The signal from a microwave source is modulated by a pulse generator through a transistor-transistor 
logic (TTLjswitch to obtain a pulsed microwave signal, which is sent into the device, and then a retrieval pulse can 
be measured after the pre-programmed retrieval time. 

(3) The third type of measurement is the coherence measurement as shown in Fig. S6:. The scheme is similar to 
the time trace measurement but with an interferometer added. The input signal is split into two branches, one of 
which is used as the reference to interfere with the output signal. The setup is very similar to the Mach-Zehnder 
interferometer commonly used in optical measurements. By varying the phase of the reference, the amplitude of the 
interference signal changes, and their relation is measured to characterize the retrieval pulse coherence. 


H. Experimental Imperfections 

The operation of the MGM requires the hybrid magnon-photon modes evenly distributed in the frequency domain, 
which will give a perfect constructive interference at the pre-programmed retrieval time T, and best suppression of 
the fringes at other times. Imperfection in the frequency distribution will result in deteriorated signal re-construction. 
This is clearly illustrated by the comparison given in Fig.[^ For a uniformly distributed spectrum, the retrieval pulse 
is very clean; while for the non-uniformly distributed spectrum, the retrieval pulse is severely distorted. 


In our experiments, we can precisely tune the small coil to obtain a near perfect gradient for the magnon modes. 
The small experimental imperfections only have very slight effects on the performance of the MGM. To investigate 
the influence of the experimental imperfection on the efficiency of the MGM, we simulated the efficiencies for various 
parameters, with random perturbations to the ideal case, where 


9j = 50 X (1 + ij), 

+ (j -^- jAidrn + AcJ^ X 


(552) 

(553) 


Here, tj ^ [—0.1,0.!] is a uniformly distributed random variable, which means the variations of magnon frequency 
and coupling strengths are within a range of ±10%. 

The results for 500 sets of different parameters are shown in Fig.[^ with the mean value of C being about 0.30, while 
the standard deviation is 0.009. For such a high imperfection up to 10%, the MGM still shows a very good ability for 
signal re-construction. Therefore, the performance of the MGM is very robust against experimental imperfections. 
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FIG. S3: a, Calculated dynamics of the MGMs with different YIG sphere numbers (N). b, Comparison of the numerical and 
analytical solutions at various coupling strengths, c & d, dynamics of the MGM with various coupling strengths and magnon 
dissipation rates, respectively. For each plot, the solid and dashed curves represent the numerical and analytical solutions, 
respectively. All parameters (beside the varying parameters) are taken from the experiments: Ns = 8, = 27r x 10 MHz, 

Au = 27 r X 10 MHz, Ka,o = 27 r x 3 MHz, Krn = x 0.72 MHz, Ka,i = i^a,o + (critical coupling), and the input pulse 

duration is tp = 20 ns. 




FIG. S4: The calculated reflection signal dynamics at (a): various input microwave detunings (ct;in —uJa)/Aoj, and (b): various 
magnon frequency gradients Act;. In the calculation, we assume all the magnon modes are identical, have uniform coupling 
strength, and are evenly distributed. Other parameters are from our experiment: Ns = 8, = 27rx 10 MHz, Au = 27rx 10 MHz, 

Ka,o = 27r X 3 MHz, Krn — X 0.72 MHz, = /^a,o + (critical coupling) and the input pulse duration is tp — 20 ns. 
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FIG. S5: a, Device schematic of the MGM consisting of eight YIG spheres, b, The intensity and direction of magnetic field 
distribution of the cavity TEno mode. Only the xy cross-section is shown because the magnetic field of this mode is almost 
uniform along the direction, c, The cosine function-like magnetic field distribution along the x direction on the cavity wall 
where the YIG spheres are located. 
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FIG. S6: Schematic of the measurement setup for the frequency spectrum of the cavity reflection (a), the time trace after a 
pulse excitation (b), and the coherence of the retrieved pulse (c), respectively. VNA: vector network analyzer; DUT: device- 
under-test; Gir.: circulator; Src.: microwave source; ED: envelop detector; PG: pulse generator; Osc.: oscilloscope; PS: phase 
shifter. 
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EIG. S7: Gomparison of the frequency spectrum and time trace for the uniformly (a) and randomly (b) distributed magnon 
mode scenarios, respectively. In each panel, the left side is the frequency spectrum and the right side is the time trace; the 
bottom rows are the results at different bias magnetic fields, while the top rows are the specific cases indicated by the dashed 
lines in the bottom rows. 
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FIG. S8: The statistics of the numerically solved memory efficiency for randomly varied parameters gj and Uj. 


















